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ABSTRACT 



We present a set of low resolution empirical SED templates for AGNs 
and galaxies in the wavelength range from 0.03 to 30/im based on the multi- 
wavelength photometric observations of the NOAO Deep- Wide Field Survey 
Bootes field and the spectroscopic observations of the AGN and Galaxy Evo- 
lution Survey. Our training sample is comprised of 14448 galaxies in the redshift 
range < 2; < 1 and 5347 likely AGNs in the range < z < 5.58. The ga laxy 
templates correspond to the SED templates presented by lAssef et al.l (120081 ) ex- 
tended into the UV and mid-IR by the addition of FUV and NUV GALEX and 
MIPS 24/im data for the field. We use our templates to determine photometric 
redshifts for galaxies and AGNs. While they are relatively accurate for galaxies 
iozjiX + -2-) = 0.04, with 5% outlier rejection), their accuracies for AGNs are 
a strong function of the luminosity ratio between the AGN and galaxy compo- 
nents. Somewhat surprisingly, the relative luminosities of the AGN and its host 
are well determined even when the photometric redshift is significantly in error. 
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We also use our template s to study the mid-IR AGN selection cri teria devel- 



oped b y IStern et al.l ( 120051 ) and iLacy et al.l (120041 ) . We find that the IStern et al. 



(120051 ) criteria suffers from significant incompleteness when there is a strong host 
galaxy component and at 2; ~ 4.5, when the broad Ha emission line is redshifted 
into the [3.6] band, b ut that it is little contaminated by low and intermediate 
redshift galaxies. The iLacy et al.l ( 120041 ) criterion is not affected by incomplete- 
ness at 2; ~ 4.5 and is somewhat less affected by strong galaxy host components, 
but is heavily contaminated by low redshift star forming galaxies. Finally, we 
use our templates to predict the color-color distribution of sources in the up- 
coming WISE mission and define a color criterion to select AGNs analogous to 
those developed for IRAC photometry. We estimate that in between 640, 000 
and 1, 700, 000 AGNs will be identified by these criteria, but without additional 
information, WISE-selected quasars will have serious completeness problems for 
z > 3.4. 



Subject headings: galaxies: active — galaxies: distances and redshifts — galaxies: 
photometry — quasars: general 



1. Introduction 



Upco ming large ph otomet ric surveys hke the ground based LSST dTvsonI 120021 ) . Pan- 
STARRS JKaisei]l2004l) . DBS Jxhe Dark Energv Survev Collaboration! I2OO5I) and UKIDSS 
( iLawrence et al.ll2007l ) projects, or the space-based WISE (IMainzer et al.ll2005l ) mission, will 
increase the number of cataloged sources in the universe by a large factor. Most objects 
always lie at the survey magnitude limit, which means that spectroscopic follow-up will only 
be possible for a small fraction of objects, and much of the science will need to rely solely on 
multi-wavelength photometric observations. One of the main tools for analyzing the multi- 
wavelength photometry is the modeling of the spectral energy distributions (SED) of the 
sources. SEDs can be used to estimate photometric redshifts and K-corrections, to study 
luminosity functions and the properties of different populations of objects, and to select 
smaller subsamples for follow-up studies. 

Galaxy SEDs have been studied using many different techniques. Theoretical models 

based on population synthe sis of stellar spectra (e.g. Bruzual fc Charlotll2003 : Fioc fc Rocca-Volmerange 
I997I : iLeitherer et al.lll999l ) have been successfully applied in many studies of galaxies. Their 
strength is that they allow the estimation of physically relevant parameters, like star forma- 
tion rates, ages, masses and metallicities solely by comparing the SED models to the pho- 
tometric observations. However, there are unknowns in the spectral properties of stars and 
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stella r populations, such as the p resence of thermally pulsat ing AGB stars (e.g. IConroy et al. 
20091 ) and the effects of binaries (lEldridge fc Stanwayll2009l ). as well as non-stellar sources of 
emission, such as dust in the interstellar medium (ISM) or active nuclei, and absorption, both 



broadly distributed and associated with stars, that alter or modi: 
issue s can be avoided by the use of empirical SEP templates (e. g. 



y the SEDs. Some of these 



Coleman. Wu fc Weedman 



1980l : lKinnev et al.lll996l : IPevriendt. Guiderdoni k Sadatlll999l : lAssef et al.ll2008[ ). However, 
most empirical templates lack the broad wavelength coverage of the theoretical ones. 

Almost all AGN SED templates are emp irical, as the physi cal processes governing AGN 
emission are not quantitatively understood. lElvis et al.l (Il994[ ) constructed the mean SED 
of Type 1 AGNs by combining spectroscopic redshifts and multi-wavelength photometric 
observations ran ging from the r a dio to the X-rays of 47 known quasars. A similar approach 
was followed by iRichards et al.l (|2006[ ) over a similar wavelength range based on 259 SDSS 
Type 1 AGNs. Both studies, however, include a limited number of AGNs, do not consider 
reddening and do not self-consistently model contamination from the host galaxy. 

In a previous paper jAssef et al.l boosl . hereafter IaOsI ) we presented a set of empirical 
galaxy templates derived fror n the extensive photo metric observations of the NOAO Deep 
Wide-Field Survey (NDWFS; IJannuzi fc Deylll999l ) Bootes field and their s pectroscopic fol- 
low up observations by the AGN and Galaxy Evolution Survey (AGES; iKochanek et al. 



m prep 



). These templates form a linear non-negative basis for the color-color space spanned 
by galaxies and accurately reproduce SEDs in the wavelength range from 0.2 to lOyum. Here, 
we extend that work in two different directions. Fi rst, we exteiid the wavelength range 
from 0.03 to 30 urn b y incorporating GALEX UV jMartin et all boosh and MIPS 24/xm 



(IWeedman et al. 



20061 ) observations, and secondly, we add a template for AGN emission 



derived in the same manner as the galaxy templates. By simultaneously deriving the AGN 
and galaxy SEDs, we self-consistently account for the host contamination. We also correct 
for the reddening of the different AGNs in our sample. Because AGES selects AGN using a 
broad range of selection methods, our AGN template is not exclusively applicable to unred- 
dened Type 1 quasars. Ea rlier versions of this extended set of templates have already been 
used bv lAtlee et al.l (120091) a s a par t of a study of the evolution of the UV upturn in elliptical 
galaxies, and by IRoss et al.l (120091 ) to study the resolved host of the z = 1.7 gravitationally 
lensed quasar SDSS J1004-t-4112. The latter study helped to verify that our models could 
correctly separate host and AGN contributions for unresolved systems. 

This paper is structured as follows. In ^we describe our photometric and spectroscopic 
data set. The algorithms we use to obtain the SED templates are discussed in §31 as well as 
the algorithm to use them for estimating photometric redshifts. In §1] we discuss our SED 
templates and their accuracy for estimating photometric redshifts of galaxies and AGNs, and 
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we apply them to study the problems of mid-IR QSO sele ction in the IRAC bands an d in the 
upcoming WISE satellite mission. In a subsequent paper ( lAssef et al.ll2009l . iPaper III ) we use 
these templates to study the IRAC-selected QSO luminosity function and its evolution with 
redshift. Throughout the text we assume a ACDM cosmology with Hq = 73 km s~^ Mpc~^, 
Qu = 0.3 and = 0.7. 



Data 



The data used in this work are an expansion over that of lAOSl . It is based on the extensive 
multi - wavelength imaging of the NOAO Deep Wide-Field Survey (NDWFS; iJannuzi fc Dey 
19991 ) Bootes field and the spectr oscopic observations of t he same field by the AGN and 
Galaxy Evolution Survey (AGES: iKochanek et al.llin prep.l ). NDWFS is a deep optical and 



near-infrared survey that covers two 9.3 deg^ regions of the sky, the Bootes and Cetus fields. 
We focus on the Bootes field, which was imaged in (3500-4750 A, peak at ~ 4000 A), R, 
I and K pass-bands to depths of approximately 26.5, 26, 25.5 and 21 AB magnitudes (5cr, 
2" diameter aperture). Extensive follow-up imaging of the field at different wavelengths has 
been obtained by several other surveys. In particular we inclu de J and Kg data fror n the 
Flamingos Extragalactic Surve y (FLAMEX: lElston et al.ll2006f) . z' data from zBootes (jCool 



20071 ) ■ MIPS 24/im data from IWeedman et al 



the GALEX GR5 release (IMorrissey et al.ll2007l). and the mid -IR IRAC observations of the 



(120061 ) . NUV and FUV observations from 



Spitzer Deep- Wide Fiel d Survey fSDWFS: lAshbv et al.ll2009l ). The IRAC Shallow Survey 



(lEisenhardt et al.ll2004l ). which we used in lAOSl . is a subset of the SDWFS data. While in 
AOSI we referred to the IRAC channels by CI, C2, C3 and C4 for the 3.6, 4.5, 5.8, and 



Sfim bands, here we will use the more common nomenclature with the channel's wavelength 
in square 



(FIRS T: iBecker et a 



NVSS: Condon et al 



1995 



Reneelink et al. 


1997 


; WSRT: c 


e Vries et al. 


2002; 


( Chandra XBootes: 


Murrav et al. 


2005) observations 



of the field that we do not directly include in our analysis. 

AGES is a redshi ft survey in the NDW FS Bootes field using the 300 fiber robotic 
Hectospec instrument (IFabricant et al.ll2005l ) on the 6.5m MMT telescope. AGES obtained 
spectra for ^ 26000 objects in the wavelength range from 3200 to 9200A. Spectroscopic 
redshifts have been measured for about 18000 galaxies in the redshift range from to 1 with 
a magnitude limit of / = 20, and about 7000 likely AGNs in the range from to 6 up to 
/ = 22.5. Because of the different magnitude limits for the samples, AGN candidates were 
pre-selected from the photometric observations. A source was targeted as an AGN if it met 
any of the following criteria: 



- 5 - 



X-Ray: The object was detected in X-rays by XBootes. The contaminants in the X-ray 
sample are primarily active st ars, along with a s mall number of very low redshift early 
and late type galaxies (see e.g. IStern et al.ll2002l for a description of the origin of X-ray 
emission in galaxies). 

Radio: The object has an associated unresolved 5a WSRT detection. Radio emission 
is generally associated with AGNs, but can also be du e to star formation at low radio 
powers (e.g. iHickox et al.ll2009t iKauffmann et al.ll2008l ). Strongly star forming galaxies 
at low redshift can contaminate this sample. 

MIPS: The object was detected in the 24/im MIPS images with a flux F > 0.3 mJy, 
appears as a point source in the optical data and has an J-band magnitude fainter than 
18 — 2.5 log(F/mJy). This J-band magnitude limit eliminates normal stars. The MIPS 
sample is contaminated by star forming galaxies that are classified as point sources in 
the optical photometry. 

IRAC: The object has IRAC colors that suggest the presence of an AGN. The exact 
color-color boundaries and magnitude limits va ry between optical point and extended 
sou rces that are discu ssed in detail in iPaper 111 , but are based on modified versions of 
the IStern et al.l (120051 ) criterion. 



The final AGES sample includes spectra for 3047, 789, 2025 and 3844 objects targeted 
as AGN by the X-ray, radio, MIPS and IRAC criteria respectively. 

After the observations, spectra were analyzed using a modified version of the SDSS spec- 
tral pipeline. Based on the template fits to the spectra, objects that showed clear signatures 
of nuclear activity were classified as spectroscopic AGNs. This classification has limits, as 
obscured AGN may lack optical signatures of activity, and the pipeline classifications em- 
phasize the presence of broad lines over the more subtle classification pro blems for narrow 



lines. For narrow line sources we also used the emission line classifications of Moustakas et al. 



(120091 ). We view these various AGN classifications as complimentary rather than focusing 
only on one. 

When possible, we have added upper limits to all objects observed, but not detected, 
in the UV, optical, near-IR and 24/im. The catalog from which we build the SED templates 
consists of all objects that have been detected, or have upper limits, in at least 8 of the 
14 bands available (FUV, NUV, Bw, R, I, z', J, Ks, K, [3.6], [4.5], [5.8], [8.0] and 24/im) 
and which have a reliable redshift measured by AGES. We eliminate objects close to bright 
stars. The minimum number of bands is required to guarantee that all SED fits have at least 
2 degrees of freedom (see Our final data set comprises 19795 objects, of which 14448 
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are nominally "pure" galaxies with no signs of nuclear activity and 5347 show some type of 
AGN signature. On average, each object has been detected in 10 of the 14 available bands, 
and Figure [T] shows the rest-frame wavelength coverage of our final sample. 

Figure |2] shows a Venn diagram representation of the 4029 photometrically selected 
AGNs in our sample. Each geometrical shape represents a different photometric selection 
criteria, and the intersections between them show the number of objects targeted simulta- 
neously by each of those criteria. The rest o f the AGNs in ou r sample are spectroscopically 
confirmed by their line ratio classification of iMoustakas et al.l (120091 1242 objects in total of 
which 1034 do not meet any photometric targeting criteria) or by the analysis of the spec- 
troscopic pipeline (2418 objects in total of which 268 were not tar geted by their photometry 
and were not classified as Type 2 AGNs by IMoustakas et al.ll2009l ). The remaining 16 active 
objects were targeted by their optical colors using a selection technique of minor importance 
in the AGES survey. Figure |2] also shows the Venn diagram representation of all 2; > 1 
objects, in order to show how the photometric selection criteria overlap in the absence of 
contaminants. Notice, however, that each selection criteria has different redshift-dependent 
biases, so the differences between the two Venn diagrams are not solely due to contaminating 
inactive galaxies. 



3. Methods 



The methods we use to deri ve th e SED templates and estimate photometric redshifts 
follow closely those described in |A08| . In this section we give a brief summary of how our 
methods wor k and we detail the changes we have made to the original algorithms. We refer 
the reader to A08 for a detailed discussion of the latter. 



3.1. Templates 

We assume that the spectrum of any object in our sample can be modeled as a linear 
combination of a small set of unknown spectral templates. We construct these templates by 
using the data described in §2J For "pure" galaxies, we assume the majority of them have 
SEDs that can be described as a linear combination of three templates: one similar to an 
elliptical galaxy (an old stellar population), one similar to a spiral galaxy (a continuously 
star forming population), and a third similar to an irregular galaxy (a starburst population). 



For objects with active nuclei, a population that was not included in the analysis of |A08 



we 



assume that every AGN SED can be described by the same spectral template with varying 
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amounts of reddening and absorption by the intergalactic medium (IGM), combined with 
the galaxy templates to describe the host. We assume the reddening follows an SMC-like 
extinction curve for A < 3300A and a Galactic extinction curve at longer wavelengths, 
assuming Ry = 3.1 for both. This is based on the observe d absence of the Galactic 2175A 
feature of Galactic extinction curves in QSO spectra f e.g. lYork et al. Ibooeh . For the SMC 
reddening law we assume the functional form derived by iGordon fc Claytoru (119981) using the 



star A zV 18, while for the Galactic reddening we assume the functional form of lCardelli et al. 



(119891 ). The IGM absorption is applied to all templates and not solely the AGN one, but it 
is most i mportant for the AGN due to their high redshifts. We model the IGM absorp tion 
following iFan et al.l ( 120061 ) for Lya and Ly/3 absorption, and lStengler-Larrea et al.l ( 119951 ) for 
Lyman limit systems. We leave the strength of absorption as a free parameter of the model 
in order to allow for variations in absorption between objects. We only include the IGM 
absorption component for redshifts above 0.8, as below this redshift Lya falls outside of our 
wavelength coverage. To force every combination of these four templates to be physically 
plausible, we require the SED of every object to be a non-negative linear combination of 
these templates. 



Our algorithm is iterative (see discussion below and §3.1 of lA08l ) and hence we must 
start from an initial guess and sequentially improve it to best fit the data at each iter- 
ation. For the galaxy templates, we start from the I Coleman. Wu fc WeedmanI (jl980l ) E, 
Sbc and Im template for th e "elliptical", "spiral" and "irregular" components, extended 
into the UV and IR bv the iBruzual fc CharlotI (120031 ) stellar synthesis models. Because 
the models of iBruzual fc CharlotI (120031 ) do not consider dust or PAH emission, which are 
critical in the mid-infrared for the star forming templates, we added ad-hoc linear combi- 
nations of the mid-infrared part of the VCC1003 (NGC 4429) and M82 SEDs obtained by 
Devriendt. Guiderdoni fc SadatI (119991 ) to the Sbc and Im templates. For the AGN tem- 
plate we consider two different starting points: an ad-hoc template constructed from dif- 
ferent observat i onal a nd theoretical considerations, and the average Type 1 quasar SED of 
Richards et al.l (120061 ). The ad-hoc template co nsists of a combinatio n of modified power- 
laws that resemble the large scale behavior of the lRichards et al.l (120061 ) Type 1 quasar mean 
energy distribution in the IR and optical, modified in the UV to fall as predicted by a simple 
accretion disk model, and with the addition of the broad-spectral lines of an AGN. We choose 
to make our in itial template a s blue in the optical as our bluest AGNs, which is significantly 
bluer than the lRichards et al.l (120061 ) SED. The shape of this template, along with the initial 
guess for our galaxy templates, is shown in Figure |3l Our results refer to this starting model 
except where noted otherwise. While the two r nodel templates produ ce very similar results, 
we prefer the ad-hoc starting point because the iRichards et al.l (120061 ) SED has considerable 
amounts of host contamination (see §4.ip . 
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Each template is divided into 300 uniformly logarithmically spaced wavelength bins 
extending from 0.03 to 30/im, hence the "low resolution". Each bin is then directly fit to 
the data using minimization over the entire calibration sample. Because this is still a 
significantly higher resolution than the width of the broad-band filters, we add a smoothing 
term (H) to the minimization that penalizes departures fr om th e smoothness of the initial 



guesses. Note that this resolution is higher than that used in |A08| . In summary, we optimize 
the function 

G = + -H. (1) 

The term is given by 

(171 TTimod \ 2 

where indices i and b enumerate objects and bands respectively, F is the observed flux, _p™°<i 
is the template model flux and ai^b is the photometric error. We defined the smoothing term 
if to be 

H = > :(log#^-log^^) , (3) 



Qk,Un Q 



where is the flux per unit frequency of SED template k in wavelength bin n (frequency 
and Qk,y„ is the corresponding value in our in itial template guess. For a detailed 



discussion of these terr as, we refer the reader to lAOSl . For galaxies, we use t] = 0.002, which 



is half of that used in IA08I due to the increased resolution, but the results are not very 
sensitive to this choice. For the AGN template we increased the value to r/ = 0.005 as the 
amplitude of the changes is much smaller than for galaxies. 

The equations for optimizing equation ([T]) are linear in all parameters but the AGN 
reddening and the IGM absorption. To solve the linear equations while imposing the condi- 
tion that the contribution of each template to each SED must be strictly non-negative and 
that at every wavelength the templates must be positive, we use a simple non-negative least 
square solver. 

To summarize, our algorithm proceeds as follows: (1) we estimate the contribution of 
each template to the SED of each object; (2) we estimate the zero point correction to each 
photometric band; (3) we sequentially optimize the templates; (4) return to step (1). Note 
that when optimizing a given template we only use objects with at least a 20% contribution 
from the template to the integrated luminosity of its best fit SED in the previous step. Since 
the 24/xm data are effectively decoupled from the rest of the dat a (see Fig. [1]), we cannot fit 



a zero-point correction to this channel. For this reason, unlike in |A08| . we do not attempt to 



remove the large scale behavior of the zero point corrections. Instead, we include a Gaussian 
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prior to penalize unnecessary modifications to the zero points, with a dispersion of 1% for 
all bands except the MIPS 24/im channel. 

We first derive the best fit galaxy SED templates without the AGN template by using 
only objects that show no signatures of nuclear activity. To be conservative, we eliminate 
all objects photometrically targeted as AGNs by AGES, all objects for which the SDSS 
pipeline finds an AGN componen t in the spectrum, and all objects classified as narrow line 



AGNs bv lMoustakas et all (120091 ) (see §2]). We then fit the AGN template by using our full 
sample while keeping the galaxy templates fixed. When determining the galaxy templates, 
we eliminate the 3% of the objects with the poorest fits, while for determining the AGN 
template we eliminate the worst 5%. These are mostly objects with either bad data points 
or galaxies with unfiagged nuclear activity. 



3.2. Photometric Redshifts 



For photometric redshifts we follow the same basic steps as those in lAOSl . We perform 



a minimization similar to that outlined in §3.11 solving all linear equations with the 
NNLS algorithm combined with a grid of reddening values for the AGN template as a 
function of redshift. The main difference is that we fix the IGM absorption to a typical 
value rather than letting it vary (see §4.21 for a discussion of this point). For optimizing the 
photom etric redshifts, w e also include a luminosity prior based on the Las Campanas Redshift 
Survey (ILin et al.lll996l ) luminosity function that only affects the galaxy templates (with and 
without and AGN component on the SED). T he sh ape and specific parameters of the prior 
probability function are discussed in detail in lAOSl . Our publicly available algorithm^!! for 
estimating photometric redshifts have been updated to refiect these changes and include the 
new templates. 



4. Results 

In this section we describe the SED templates resulting from applying the algorithm 
discussed in §3.11 to the data set described in §2] and we also discuss their application to 
determining photometric redshifts for AGNs and galaxies, and to study the mid-IR selection 
of QSOs in the IRAC bands and in the upcoming WISE satellite mission. 



Codes and templates are available at 



www . astronomy . ohio-state . edu/~r j assef /Irt 
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4.1. Templates 



The final templates are shown in Table [T] and in Figure [31 Table [2] shows the absolute 
magnitudes of the templates as a function of redshift in various photometric bands including 
those of our sample. Figure |3] also shows the initial guesses from which the process star ted. 
The galaxy templates show significant departures from the initial ones, just as in lAOSl In 



particular, the elliptical template is significantly redder in the UV and somewhat bluer in 
the mid-IR. The spiral-like template has decreased UV and a increased dust/PAH emission 
relative to the 1.6/im stellar emission peak. The Im template shows a somewhat decreased 
UV and mid-IR continuum relative to the near-IR peak. The departures from the initial 
templates are caused by a combination of changes in the balance of the underlying stellar 
populations and by their inability to reproduce the SEDs of galaxies, either because of 
missing physics or, in the mid-IR, because of the ad hoc addition of dust emission used to 
create them (see §3.11 for details). 

The AGN template based on our ad-hoc model also shows significant differences from 
the initial guess. In particular, the optical is bluer and the mid-IR is redder. At longer 
wavelengths (> 15/im), the template shows a steep rise, but this is likely to be an artifact 
caused by the lack of z = quasars to constrain this wavelength region. At wavelengths 
shortward of Lya, the templates cannot be constrained, as the observed SEDs are completely 
dominated by IGM absorption. Since the flux in X-rays is generally signif icantly lower than 



the UV emission in AGNs, as determined by the values of aox (e.g. iTananbaum et al. 
19791 : IStrateva et al.l l2005l ). we know the spectrum must turn down at some wavelength 
shortward of Lya. For this reason, we choose in the rest of the paper to estimate "bolometric" 
luminosities by integrating our AGN template only at wavelengths longward of Lya. 



Using the average SED of Type 1 quasars derived by [Richards et al.l (120061 ) as a starting 
point for our AGN template yields similar results, albeit with some interesting differences. 
The AG N template resulting f rom this starting point is shown in Figure HI as compared to the 
original [Richards et al.l ( l2006l ) template and our "standard" template. Relative to the initial 
template, the converged template is somewhat brighter in the UV relative to the optical and 
near-IR wavelengths. The mean color is a combinati on of extinction and any dependence of 
the SED on black hole mass and accretion rate (e.g. lYip et al.ll2004l ). While our algorithm 



corrects for the former, it can mimic the latter using small amounts of reddening. In the 
optical wavelengths, the contribution from broad emission lines to the iRichards et al.l (|2006[ ) 
template is heavily smoothed by the averaging of the broad band photometry used to build 
it, however the output template to our algorithm has developed much sharper emission line 
features. The minimum in the near-IR of the converged template has shifted to a shorter 
wavelength, as also seen in our "standard" template. This is likely due to our self-consistent 
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treatment of the host. At ~ 4.5/im there is a depression in the SED. It is unclear if this is 
a real feature or not, as it does not appear in our standard AGN template. The differences 
at the longest wavelengths and wavelengths shorter than Lya are not significant due to the 
lack of data needed to constrain these regions (see Fig. [T] and the discussion in the previous 
paragraph). When combined with the host gala xy templates, both AG N templates work well. 



suggesting that the main difference is that the [Richards et al.l (120061 ) template intrinsically 
has more host contamination than our standard template, which can be compensated for by 
changing the amount of host galaxy template added to fit an object. This can be seen in 
Figure [5|, which shows th a t the values of LAGN/-Z^Host are, in average, ~ 35% larger for the 
converged Richards et al. AGN template than for the standard template. 



In Figure Owe show the best fit extinctions to the AGNs in our sample for our standard 
AGN template. Objects have been divided according to the criteria by which they were se- 
lected and objects targeted by multiple criteria are shown in all of the relevant panels. In all 
cases, there seems to be a maximum observed reddening that increases with increasing AGN 
luminosity relative to their hosts. This can be understood as a combination of three effects. 
One is that a more heavily reddened AGN needs to be more luminous relative to the host in 
order to meet many of the AGES AGN targeting criteria. The second effect is the inability 
of our templates to detect highly reddened AGNs with bright hosts and, if detected, to de- 
termine the reddening parameter with precision. If the host dominates the UV/optical/NIR 
wavelengths, reddening values above a certain threshold will not significantly affect the 
until they start to modify the mid-IR colors. Otherwise, our algorithm defaults to using 
the smallest possible reddening value. Finally, our sample is magnitude limited, and hence 
the more reddened, optically fainter objects can be found in a smaller physical volume than 
the less reddened, optically brighter ones. The first problem primarily affects the IRAC and 
MIPS targeting criteria, while the second and third ones affect all cases. Notice, however, 
that the X-ray selected objects are relatively unaffected by these biases. In Figure E] we 
also see that the AGNs which are brightest in comparison to their hosts show a range of 
small, but positive, extinctions. This may be due to our algorithm using small amounts of 
reddening to compensate for changes in the SED with BH accretion rates and masses, as 
discussed above. 



Figure [3] also shows the galaxy templates derived in lAOSl . The elliptical and irregular 
templates show significant changes only in the UV and mid-IR, as would be expected from the 
addition of the GALEX and 24/im MIPS data. The spiral template also shows departures 
at those wavelengths, but differences are also seen in the near- to mid-IR region. In our 
model the templates form a non-orthogonal basis that, due to the restriction to positive 
combinations, determine the boundaries of the color space occupied by galaxies. Changes 
in the templates are really changes in the boundaries of the permitted color space. We 
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expect that the differences between the new and old templates over the wavelength range 
covered by the old templates can largely be explained by making the new/old templates 
linear combinations of the old/new templates. Note that the coefficients of the "ro tation " 



is 



can be negative. The best fit combination of our new templates to the E template of I AOS 
given by the coefficients (1.02, -0.01, 0.00), to the Sbc one by (-0.51, 1.51, 0.00) and to the 
Im one by (-0.04, 0.34 and 0.70), where each coefficient represents the contribution of our 
new E, Sbc and Im tem plate s respectively. These combinations reproduce very well the full 
wavelength range of the lAOSl templates, showing that the differences with our new templates 
correspond to a simple rotation of the basis vectors rather than a fundamental modification. 

Figures [7] and [8] show typical fits to AGNs and galaxies using the "standard" AGN 
template. These objects have values close to the median of their respective samples 
(median = 1.35 for galaxies and 2.60 for AGNs) and have been selected to show a range 
of different scenarios. The bottom panel of each Figure also shows the object with the 
90**^ percentile worst of each sample. Although small discrepancies are present at some 
wavelengths, the agreement is still generally very good. The vast ma jority of the objects in 



our sample are well described by our four templates. The converged [Richards et al.l (120061 ) 
template gives similar values of median (1.30 and 2.61 for galaxies and AGNs respectively) 
to our standard template as we would hope, given that our results should depend little on 
the priors. 

In the following sections, we will explore different applications of our templates. We will 
study the accuracy of photometric redshifts from these templates for galaxies and AGNs, 
and the colors of these objects i n the mid -IR bands of both IRAC and the upcoming WISE 



mission. In a subsequent paper (jPaper 111 ) we examine the AGN luminosity function and its 
evolution with redshift. Our templates also allow for the study of host properties in AGNs 
but this will be explored in an upcoming paper. 



4.2. Photometric Redshifts 

In this section we study the application of these templates to the estimation of photo- 
metric redshifts for galaxies and AGNs. To quantify their accuracy, we derive photometric 
redshifts for all the sources in the catalog used to build the templates, each of which has 
a measured spectroscopic redshift. As mentioned previously, for photo-zs we do not allow 
the IGM absorption strength parameter to vary, as allowing variations worsens the redshift 
accuracy, and we include a luminosity prio r based on the Las Campanas Redshift Survey 
r-band luminosity function (jhin et al. 1996 ). 
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We divide our catalog into three distinct groups: "pure" galaxies, extended AGNs 
(optically resolved galaxies with AGN signatures), and point source AGNs (objects where 
the host galaxy was not optically resolved), based on the AGES targeting criteria (see 
The difference between the two AGN groups is a combination of redshift and the relative 
brightness of the host and the AGN. Photometric redshifts for all groups are estimated in 
the same way, using all four templates and using a photometric redshift range of < 2; < 6. 
To study the accuracy of the photometric redshifts, we calculate (i) the standard dispersion 
az/{l + z) defined as 



[l + z) 



1 ^ 



1 + zi 



(4) 



where Zp and z^ are the photometric and spectroscopic redshifts respectively, (ii) the median 
offset of Zp — Zs, (iii) the ranges of \zp — Zs\/ {)- + z^) encompassing 68.3%, 95.5% and 99.7% of 
the distribution, and (iv) the dispersion Az which is equivalent to ^^/(l + z) but uses only 
the 95% objects with the most accurate photometric redshift estimates in order to eliminate 
outliers. 

We note that including the MIPS 24/im data slightly worsens the photometric redshift 
accuracies for all groups. This is in contrast to including the mid-IR data, which improves the 
accuracy. Depending on the object type, including the MIPS data increases the dispersion (as 
measured by Eqn. |3]) by factors between 1.07 (galaxies) and 1.21 (point source AGNs), and 
the 95% clipped accuracy by between 1.05 (extended AGNs) and 1.09 (point source AGNs). 
This is probably symptomatic of a scatter in the observed SEDs larger than permitted by 
our templates, even when they fit the mean sample characteristics well. This may be caused 
by a small number of LIRGs and ULIRGs that should be modeled as a different population 
at these wavelengths. For this reason we exclude the MIPS photometry when estimating 
photometric redshifts for the rest of the paper. 

The results of our photometric redshift estimates are summarized in Table [31 Figure 
|9] shows a comparison of Zp against Zs for the galaxy sample, while Figure [10] shows it for 
the extended and point source AGNs. For galaxies, the contours closely follow the Zp = Zg 
line, showing that our photometric redshifts estimates are accurate and that most galaxies 
can indeed be described by our 4 template model. The 95% clippe d redshift accuracy is 
= 0.041. These results are compatible with those presented in lAOSl . although the accuracy 
is slightly worse. This is largely due to expanding the redshift search range from 0-1 to 0-6 
in order to encompass the AGN population, even when it i s not physical to do so . It is worth 
noting that this accuracy is comparable to that found by iBrodwin et al.l ( 120061 ) for galaxies 
in the NDWFS Bootes field using a hybri d template fitti i ig an d neural network approach 
and an earlier version of the AGES data. iBrodwin et al.l ( 12006[ ) found (^^/(l + z) = 0.105 
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and Az = 0.047, about 20% better and 15% worse, respectively, than the accuracy we find. 

The photometric redshift accuracies for AGNs are much worse. For extended AGNs, 
the dispersion is 20% larger than for galaxies (Az = 0.050) and for point source AGNs it 
is a factor of ~ 4.5 w o rse (A z = 0.18). This lower accuracy is similar to the results of 



Rowan-Robinson et al.l (120081 ) for QSOs in the SWIRE survey. The problem can be largely 
attributed to the lack of prominent broad-band features in the SED of AGNs as compared 
to galaxies. For many objects, the AGN continuum combined with that of the host produces 
an extremely fiat SED, and for these objects the photometric redshift estimate is dominated 
by secondary factors like photometry errors and the galaxy luminosity prior we apply to 
the hosts. The latter has the effect of assigning severely underestimated redshift values to 
many sources, although removing it improves the median shift but worsens the accuracy. 
Figures [TT] and [12] show the average best fit SEDs for the point source AGNs with "good" 
and "bad" photometric redshifts (defined by \zp — Zs\ < 0.2 and \zp — Zs\ > 0.5 respectively) 
as a function of the estimated reddening of the AGN. It is very clear that the objects with 
good photometric redshifts have significantly stronger host components than those with bad 
estimates. In particular, the objects with good redshifts appear to have a distinctive infiexion 
point close to 4000 A which probably constrains the photometric redshift. 

To further explore this point, we show in Figure [13] the difference between the photomet- 
ric and spectroscopic redshifts as a function of the ratio between the bolometric luminosities 
of the AGN and host components for all objects classified as point source AGNs. Note that 
we have not included the galaxy luminosity prior when fitting the SED at Zp (although Zp is 
calculated with the prior), as it is not included for the fits at Zs to which we are comparing. 
Independent of whether we estimate the ratio at Zg or Zp, the photometric redshift accuracy 
is markedly better when the best fit AGN component is fainter than the host component. 
For Lagn < -^Host, the photometric redshift dispersion is Az = 0.041 for extended AGNs and 
Az = 0.081 for the point source ones, while when Lagn > -^^Host, the dispersions rise to 0.145 
and 0.224 respectively. Because there is a marked dependence of the photometric redshift 
accuracy on the host galaxy strength relative to the AGN SED, an IRAC color dependence 
is also expected. Figure [H] shows the IRAC color-color diagram where the points have been 
color-coded according to their photometric redshift accuracy. Most of the objects with bad 
redshift estimates lie at colors redder than [3.6]-[4.5] = 0.6. 

From Figure [13] it is possible to see that the LAGN/-^^Host ratio is little changed whether 
we use the true or photometric redshift of the source. This is more clearly shown in Figure 
[TS]and has two interesting implications. First, since this ratio is not highly dependent on the 
photometric redshift, it can be used as a sign of the photo-z's likely reliability even in the 
absence of spectroscopic confirmation. And second, the ability of our templates to measure 
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Lagn / -Z^Host is not dependent on whether the true redshift of the source is known, and hence 
can be used to target (essentially) high Eddington ratio AGNs using purely photometric 
redshifts. 

At the highest redshifts of our sample {z > 4), photometric redshifts are affected by a 
color degeneracy between the Lyman break at high redshift and the Balmer break at low 
z, mainly caused by their lack of detection in the GALEX and Bw bands. It is clear from 
Figure [TD] that if this degeneracy is broken, by, for example, a combination of deeper UV 
imaging and stronger upper bounds, the photometric redshift accuracy for these objects can 
be very good. If we recalculate the photometric redshifts for these objects but only allowing 
for Zp > 1, we obtain an accuracy of (Jz/{1 + z) = 0.072. Note that this degeneracy is 
not related to the fact that we fixed the IGM absorption, as the accuracies for these high 
redshifts worsens if we allow the IGM absorption strength to vary. 

We have shown the lack of accuracy in the broad-band photometric redshifts of AGNs 
is a product of the lack of spectral features in their SEDs. Regardless, it is important to 
understand if these photometric redshifts really just lack precision or truly lack accuracy. 
Based on the probability distribution function (PDF) of each object, we have calculated 
68.3% and 95.4% errorbars on the photometric redshifts. The error bars are asymmetric 
and have been defined to include the given percentage of the PDF to either side of the 
maximum. For inactive galaxies we find that the average difference between the la {2a) 
upper and lower limits of the photometric redshifts is 0.08 (0.16), while we find 0.10 (0.21) 
for optically extended AGNs and 0.25 (0.54) for optical point source AGNs. This shows 
that our photometric redshifts for AGNs lack the precision of those of inactive galaxies and 
extended AGNs. If we recalculate az/{l + z) (Az) limited to objects where the 68.3% range 
of the photometric redshifts is less than 0.1, the range of the optically extended AGNs, we 
find 0.062 (0.036) for inactive galaxies, 0.109 (0.036) for optically extended AGNs and 0.253 
(0.118) for optical point source AGNs. The lack of improvement in these statistics for the 
optical point source AGNs suggests that the lack of spectral features in the SED of the AGNs 
not only worsens the precision, but also the accuracy of the photometric redshifts. 

The low accuracies of AGN photometric redshifts we have obtained at z < 4 can be gen- 
erally extrapolated to all SED fitting algorithms relying solely on broad-band data in the UV 
to mid-IR wavelength regime, as it is caused by the lack of strong spectral features inherent 
to the SEDs of AGNs and not to a method in particular. Nonetheless, other approaches that 
do not rely solely on broad-band SED fits have proven to be significantly more successful. 
In particular, the addition of narrow and medium photom etric bands observ ations can help 



constrain the photometric redshifts for AGNs, as shown by lWolf et al.l (|2003[). who o b taine d 



a dispersion of a = 0.03 for AGNs in the COMBO-17 survey, and bv ISalvato et all fl2009f ). 
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who obtained a /{I + z) ~ 0.015 for QSOs with 1*^^ < 24.5 in the COSMOS field. Narrower 
bands significantly enhance the photometric effects of emission lines. Estimates based on 
neural networks can also provide accurate photometric redshifts for AGNs. In particular, 
Brodwin et al.l (120061 ) used this approach to obtain redshifts for QSOs using an earlier ver- 
sion of the photometry used in this paper for the NDWFS Bootes field. We have compared 
their redshift estimates with ours for a sample of 2290 active objects with SEDs dominated 
by the AG N. While our SEP fit ting method gives o"/(l + z) = 0.18, the neural network 
approach of iBrodwin et al.l (120061 ) has an accuracy cr/(l + z) = 0.13, about 25% better. 



4.3. Mid-IR AGN Selection 



Lacv et al. 


2004; 


Stern et al. 


2005: 



naras et ai.llzuuDi ). as 

their properties at these wavelengths are typically very different from those of stars and 
galaxies. Mid-IR selection methods also avoid some of the biases against extincted AGN or 
AGN with z ~ 3 found in optical surveys. In this section we use our galaxy and AGN SED 
templates to study the selection of active galaxies by their Spitzer IRAC colors and possible 
selection biases as a function of redshift and redd ening of the centra l en gine. In par t icular , 



we will focus on the selection criteria developed by lStern et al.l (120051 ) and lLacy et al.l (120041 ) 



as they are the most commonly used AGN selection methods based on IRAC colors. 

Figure [T6l shows the mid-IR color space of the SDWFS survey (only object s with J < 21.5 
to limit the contamination by unflagged stars) along with the IStern et al.l (120051 ) quasar 
selection criterion and the color tracks of our templates (2; < 3 for the galaxy templates and 
2; < 10 for the AGN). All low redshift galaxies have very similar [3.6]-[4.5] colors, as this 
is dominated by their stellar emission. However, they populate two different regions of the 
[5.8]-[8.0] color, where the difference is dust/PAH emission in star forming galaxies and its 
absence in ellipticals. As the redshift increases, the colors converge and follow very similar 
tracks, as the bands are redshifted into the predominantly stellar emission at wavelengths 
shortward of ~ 5^m. The IRAC colors of AGNs are, on the other hand, dominated by hot 
dust emission at low redshifts/long wavelengths, and by emission of the accretion disc at 
higher redshifts/short wavelengths. As a result, AGNs populate a very distin ct region of 



the IR AC color-color space from low and intermediate redshift galaxies, and the lStern et al. 



(120051 ) or iLacy et al.l (120041 ) selection criterion approximately delineate this region. Notice 
that, as shown in Figures [16] and [T71 high redshift galaxies {z > 3) occupy a v ery similar 



region to AGNs in the IRAC colors, and contaminate both selection criteria. lYun et al. 



(120081 ) has also shown that this is the case for submillimetre-bright galaxies at z > 2. 
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However, since AGNs are commonly much more luminous than galaxies, these issues are 
only significant for faint sources. 



The IStern et al.l ( 120051 ) selection criterion does a remarkably good job of separating 
the low and intermediate redshift galaxies from the AGN. Its only significant problems are 
the "blue loop" out of the region at 2; ~ 4.5, where the Ha line is redshifted into the 
[3.6] band, and the requir ement that the AGN, rather than the host, dominates the mid-IR 
SED (jGorjian et al.l 120081 ). A survey with mid-IR completeness problems near the z ~ 4.5 
loop can be easily handled by including bluer mid-IR objects that also show evidence for 
a Lyman limit in the optical. At z ~ 7 our unreddened AGN template predicts a similar 
problem when H/3 is redshifted into the [3.6] band, and the n iid-IR colors resemble those 
z ~ 1 — 2 galaxies. Moderately reddened AGN stay inside the IStern et al.l ( 120051 ) selection 
boundaries at all redshifts, as illustrated in Figure [TBI for E{B — V^ = 0.4. With this amount 
of reddening, when Ha is redshifted into the [3.6] channel, the colors reach the selection 
boundary but are always inside of it. Note, however, that there is a host dependence on 
the amount o f redd ening necessary to stay inside the selection boundaries. As shown by 
Hickox et al.l ( l2007l ) (see also the top left panel of Fig. E]), more highly reddened objects can 
more easily exit the AGN selection boundaries due to host contamination, as it is harder to 
meet the requirement that the AGN dominates the SED. 



The mid-IR AGN selection criterion defined by iLacy et al.l (|2004| ) does not suffer from 
the z ~ 4.5 incompleteness problem. However, it does suffer from large amounts of contam- 
ination by low reds hift star forming g alaxies. Figure [T7| shows the SDWFS sources in the 
color space used by iLacy et al.l ( 120041 ) to define their AGN selection criterion. The AGN 
colors stay inside the selection boundaries throughout the whole redshift range from 0-10. 
However, star forming galaxies characterized by the Sbc template are also contained within 
the selection region at 0.25 < z < 1. Note that this also means the requirernent th at the 
AGN dominates the SED is somewhat more relaxed than for the IStern et al.l ( l2005l ) selec- 
tion criterion. A redder S's.s/S's.e limit would eliminate most star forming galaxies, but at 
the cost of introducing the bias against high redshift QSOs. Hence, the differences between 



the IStern et al.l (|2005[ ) and iLacy et al.l (|2004[ ) AGN selection criteria can be regarded as a 
trade-off between comple teness and g alaxy contamination. Figure [T7I also shows an updated 
version of this criteria by lLacy et al.l ( 120071 ). on which the vertical boundary h as been shifted 
by 0.1. This modified boundary keeps the completeness of iLacy et al.l (|2004[ ) and somewhat 
reduces the galaxy contamination. 



Recently, iRichards et al.l ( 120091 ) used IRAC observations of SDSS objects to compare 
both of the se lection crite ri a, an d their results are in very good agreement with ours. They 
find that the IStern et al.l (120051 ) selection criteria is much less contaminated by inactive 
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galaxies in comparison to that of iLacy et al.l (|2004| ). but that the lStern et al.l ( 120051 ) selection 
is biased against AGNs in the redshift range 3.5 < z < 5. We find that this problem is 
restricted to a somewhat narrower re dshift range. In an ear lier work, iDonley et al.l (|2008[ ) 
used objects in the GOODS-S field ( iGiavalisco et al.l 120041 ) to compare the reliability of 
th ese two AGN selectio n crite ria in addition to the pow er - law galaxies selecti on method 
of lAlonso-Herrero et al.l (120061) and iDonley et al.l (120071 ) . iDonley et al.l ( l2008l ) concludes 
that the Stern et al.l ( 2005) selection dia gram has more overlap with color tracks of inactive 
galaxies than that of iLacy et al.l (l2004j ). suggesting that AGN selection by the former is 
subject to more contamination. This is possibly due to including significantly fainter galaxies, 
but the dominant cause seems to be that their SED models exaggerate the color space 
populated by galaxies. As can bee seen in t heir F igure 6 for objects with 0.25 < z < 0.75, 
few non-AGN actually lie in the lStern et al.l (|2005[ ) region even though the SED tracks would 
allow them to do so. 



4.4. Colors of Galaxies and AGNs in WISE 



The Wide-Field Infrared Survey Explorer (WISE; iMainzer et al.ll2005[ ) commenced an 
all-sky survey in the mid-IR in late 2009. The survey will be made in 4 bands with effective 
wavelengths of approximately 3.4, 4.6, 12 and 22/im, estimated by assuming a source with 
flux Fi, oc z/~^. The two shortest wavelength bands are similar to the [3.6] and [4.5] channels 
of IRAC while the longest wavelength band is similar to the MIPS 24/im channel. Following 
the notation we have adopted for the IRAC bands, we will refer to the WISE bands by [3.4], 
[4.6], [12] and [22]. The average survey depth (8 repeated observations, 5a) will be 0.12, 0.16, 
0.65 and 2.6 mJy, respectively. For comparison, the average depth (6" apertures, 5a) of the 
SDWFS survey in the two shortest wavelength IRAC bands are 5.2 and 7.2 /iJy respectively 
(lAshby et al.ll2009l ). about 20 times deeper than the average WISE fleld. 



We can combine our object samples and templates to synthesize WISE observations. 
For this purpose we use the measured WISE fllter curves, including the detector response 
(Wright 2009, private communication). These are shown in Figure [H Table |2] shows the AB 
absolute magnitudes of our templates in the WISE bands as a function of redshift. Figure 
[T8] shows the WISE colors as synthesized by our templates for all SDWFS sources that will 
be brighter than the WISE [3.4] and [4.6] magnitude limits and have / < 22.5. Figure [T^ 
shows the color tracks of our A GN template (0 < ^ < 6), our star forming galaxy templates 
(0 < z < 2) and the GRASIL Jsilva et ahlEggsl ) Arp 220 SED {0 < z < 3). The latter is 
also shown in Figure [18] along with the color tracks of our E template (0 < z < 2). Note 
that the low redshift {z < 0.5) [12] and [22] fluxes of AGNs are not reliable because there 
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are too few pure AGNs in our sample at these low redshifts to stably determine the AGN 
SED at these wavelengths. The apparent lack of pure elliptical galaxies in Figure [18] is an 
artifact. Tiny contributions from the star forming SED templates significantly modify the 
colors of such galaxies in the two long est wavelength ban ds, whether these contributions are 



real (e.g. 24/xm excess red galaxies of iBrand et al.ll2009l ) or simply an artifact due to model 
uncertainties. 

Of the 3973 sources in Figures [T51 and 3694 have spectroscopic redshifts and we use 
photo-z's for the rest. We have eliminated all photo-z objects with fits having > 50 (63% 
of the photo-z sources) in order to limit contamination by stars. Since the WISE survey 
is significantly shallower than SDWFS, essentially all the AGNs at the WISE depths have 
spectroscopic redshifts and our results are little affected by any problems with their photo- 
metric redshifts. Objects targeted by AGES as AGNs are marked in Figure ITS] according to 
the criterion by which they were selected. Objects targeted by multiple criteria are marked 
only by one criterion, where the priority ordering was IRAC, MIPS, X-ray and, lastly, radio 
selection. IRAC selected AGNs clearly populate a distinct region of the WISE color-color 
diagram, well separated from galaxies. 

Non-IRAC targeted X-ray and radio candidates have a much larger scatter and tend to 
populate a region that extends between the colors of IRAC selected AGNs and normal star 
forming galaxies. This is explained by the requirement that IRAC selected AGN must have 
higher Eddington ratios, so that color is dominated by the AGN in the rest-frame mid-IR. 
Lower Eddington ratios have colors contaminated by the host, and hence have mid-IR colors 
increasingly dominated by stars and cold dust emission from star formation. X-ray sources 
are little affected by h ost contamination a. nd so are still found when the mid-IR colors are 



dominated by the host ( iGorjian et al.ll2008l ). while the WSRT radio survey is deep enough to 
detect radio emission associated with star formation as well as AGN. Note that for z > 0.5 
the host contamination needed to move Type 1 AGNs out of the selection region is not just 
a function of its SED and the Eddington ratio of AGN, but also of reddening, as both, the 
[3.4] and [4.6] channels can be affected by extinction. 



In an analogous way to the AGN selection region defined by IStern et al.l (120051 ) . we 



propose the following color-color criterion for selecting AGNs in WISE observations, 

[12] - [22] > 2.1 (5) 

[3.4] - [4.6] > 0.85 (6) 

[3.4] - [4.6] > 1.67 ([12] - [22]) - 3.41, (7) 



as shown in Figures [18] and [19] Like the IStern et al.l (l2005l ) selection method, it will have 
little contamination from low and intermediate redshift galaxies. The left boundary is chosen 
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because there do not seem to be any AGNs at redder colors either observationally or based on 
our templates. The other two boundaries are selected as a compromise between maximizing 
the number of WISE selected AGNs and limiting the contamination by intermediate redshift 
galaxies and ULIRGS, were we have characterized the latter by the SED of Arp 220. In 
particular, the rightmost boundary separates AGNs from low {z < 0.2) and high {z > 2) 
redshift ULIRGs. ULIRGs are relatively rare and interesting, so the simple and deeper 
criterion for selecting AGNs of 

[3.4] - [4.6] > 0.85 , (8) 

(the dashed line in Figure [T8|) may be all that is needed. Note that all WISE selection 
criteria have significant problems at z > 3.4, starting when Ha emissio n is redshifted into 



the [3.4] band, analogous to our earlier discussion of the lStern et al.l (120051 ) selection criterion. 
However, in this case, reddened AGN do not re-enter the selection region at higher redshifts. 
WISE data will need to be merged with other information to separate high redshift QSOs 
from lower redshift galaxies. 

Inside the full selection color region described by equations (E]), (jSD and ([7j), we find 
140 Bootes objects flagged as AGN by one of the methods described in §2] and that are 
bright enough to be detected by WISE in all of its channels. Since the full NDWFS Bootes 
field survey area is approximately 9 deg^, scaling to a full sky survey area we predict that 
WISE should find 640, 000 AGNs, corresponding to a number density of 16 deg~^, using our 
color selection criteria. If only the criterion of equation dHJ is applied, we find 375 objects 
flagged as AGN and bright enough to be detected in the two shortest wavelength channels of 
WISE. Again scaling to the full sky area of WIS E, we predict it sho uld find 1, 700, 000 AGNs 



corresponding to a number density of 42 deg~^. lAshby et al.l (120091 ) estimated a high surface 
density but used a bluer cut of [3.4]-[4.6] which will introduce significant contamination by 
ULIRGS and star forming galaxies at low and intermediate redshifts. 

These surface densities are higher than those of the SDSS survey, which finds 11.2 deg~^ 
with z < 2.3 and i < 19.1 and 1.03 deg~^ and i < 20.2, but the differences are not a simple 
change in this numbers. Figure [20] shows the z-band magnitude limit corresponding to the 
WISE band detection limits for several models of AGN SEDs. For our pure AGN template 
without a host or extinction, the effective i-band depth of the WISE sources changes strongly 
with redshift, reaching a minimum at z ~ 4.5 and then rising again. The drop at low redshift 
is due to the v-shaped structure of the SED (see Fig. [3]) - the ratio of optical to mid-IR fluxes 
rises with redshift until the minimum at 1/xm enters the mid-IR band passes at 2; ~ 4, while 
at the same time the IGM absorption begins to remove optical flux. Adding a typical level 
of host emission markedly reduces these effects by filling in the minimum of the SED, to 
make the effective optical magnitude limit of WISE for the [3.4]/[4.6] bands deeper than 
i = 19.1 at almost all redshifts. Adding a small amount of extinction to suppress the optical 
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relative to the mid-IR has similar effects. Thus, the WISE sample will complement the SDSS 
samples. 



Conclusions 



We have created an optimized basis of low resolution SED templates for AGNs and galax- 
ies in the wavelength range from 0.03-30/im derived from the extensive multi-wavelength 
observations of the NDWFS Bootes field and the AGES spectroscopic survey. This basis of 
templates consists of three galaxy SED templates that resemble elliptical, spiral and irregu- 
lar galaxies, and an AGN SED template with variable reddening. Our model also includes 
IGM absorption of variable strength for high redshift sources. The templates and source 
codes needed to synthesize other bands and determine photo-zs are publicly availabl^. 

In this paper we have investigated three applications of our templates, namely the 
estimation of photometric redshifts for galaxies and AGNs, the reliability of AGN selection 
in the IRAC bands, and the prediction of the c olors of t hese objects in the WISE survey 
as a function of redshift. In a subsequent paper (jPaper 111 ) we have used them to study the 
luminosity function of mid-IR selected quasars from z = — 5.85. 

The photometric redshift accuracy of our tem plate s for galaxies is good and comparable 



to other methods in the literature and to those in lAOSl . For AGNs, however, the accuracy is 
much lower, in particular for those that have little contribution f rom the host galaxy. This 
is in agreement with the results of iRowan- Robinson et al.l (120081 ) for QSOs in the SWIRE 
survey. The photometric redshift accuracy is highly correlated with the relative luminosities 
of the AGN and host components (i^AGN/-^Host)) in the sense that the accuracy is lower for 
higher AGN fractions. In this limit, the SEDs of AGNs are fiat, leaving no strong features for 
determining the photo-zs in broad-band photometry. When only objects where Lagn < -^Host 
are considered, the photometric redshift accuracy is comparable to that of normal galaxies. 
We find that the estimate of i^AGN/-^Host is robust even when the photo-z is inaccurate. 

We have also used our templates to study the color distribution of galaxies and AGNs 
in the IRAC bands as a function of redsh ift and AGN redd ening, and how this affects the 
selection criteria of 



that the Stern et al. 



icicv et al.l ( I2OO4J ) and IStern et al.l (120051 ). In particular, we have shown 
( 2OO5I ) criterion suffers from significant incompleteness at 2; ~ 4.5, due 



to the broad Ha line being redshifted into the [3.6] band, but that it is little contaminated 
by low and intermediate redshift galaxies. We have also shown that the selection criterion 



^Templates and codes available at 



www . astronomy . ohio- state . edu/~r j assef /Irt 
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of iLacy et al.l (120041 ) is less subject to these sources of incompleteness but it is heavily 
con taminated by low redshift star-forming galaxies. Moreover, shifting the boundaries of 
the iLacy et al.l ( 120041 ) selection region to limit the conta mination by galaxies creates similar 
incompleteness issues to those of the IStern et al.l (120051 ) criterion. The differences between 
the two selection methods can be regarded as a trade-off between completeness and galaxy 
contamination. 

We defined two simple AGN selection criteria for the WISE satellite mission by synthe- 
sizing the WISE bands for all sources in the much deeper SDWFS survey of the Bootes field 
using our template models. The absolute magnitudes of our templates in the WISE bands 
as a function of redshift are provided in Table [2l A restrictive criteria using all 4 WISE 
bands would identify 140 AGN in the Bootes field, while one based only on the 2 shorter 
wavelength bands would identify 375. The main differences are that the more restrictive 4 
filter criterion would avoid contamination by ULIRGs but is shallower because it requires 
detection in the less sensitive longer wavelength bands. Extrapolating from the ~ 9 deg^ 
Bootes field to the full sky, the two criteria would identify 640, 000 and 1, 700, 000 AGNs 
respectively but largely with z < 3.4. Additional information is needed to identify objects 
with z > 3.4. These criteria shoul d work also reasonably well in regions of high st ellar density 



but l ittle star formation (see e.g. iKozlowski &: Kochanekll2009l on applying the IStern et al. 



20051 approach to find QSOs behind the Magellanic Clouds). 
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Table 1. SED Templates 



A(^m) 


AGN (xlQ-ii) 


(er 

AGN 2 (xlO-i^) 


g/s/cm^/Hz) 
E (x 10-16) 


Sbc (xlO-i^) 


Im (xlO-i^) 


0.0303 


0.9662 


5.5387 


0.9659 


1.9502 


1.1327 


0.0311 


0.9838 


5.7997 


1.0250 


2.1322 


1.2553 


0.0318 


1.0008 


6.1105 


1.0829 


2.2252 


1.3064 


0.0325 


1.0153 


6.4724 


1.1465 


2.3982 


1.4211 


0.0333 


1.0299 


6.8141 


1.2075 


2.5162 


1.4919 



Note. — Electronic table that presents the flux per unit frequency Fi, of our templates. 
The column marked as "AGN" contains our standard AGN template while that marked 



"AGN 2" contains the SED derived by using the iRichards et al.l (120061 ) SED as a starting 
point. Templates are normalized to be at a distance of 10 pc and to have an integrated 
luminosity between the wavelength boundaries of 10^^ Lq. 



Table 2. SED Templates Absolute Magnitudes 



2 Template NUV FUV Bw B V R I z J H Ks K [3.6] [4.5] [5.8] [8.0] 24/im [3.4] [4.6] [12] [22] DM 



0.0 


1 


-15.78 


-15.69 


-15.48 


-15.36 


-15.34 


-15.69 


0.0 


2 


-12.00 


-13.63 


-18.51 


-18.59 


-19.50 


-20.13 


0.0 


3 


-14.66 


-15.86 


-18.26 


-18.28 


-19.00 


-19.59 


0.0 


4 


-17.56 


-18.05 


-19.24 


-19.22 


-19.57 


-19.88 


0.1 


1 


-15.87 


-15.76 


-15.63 


-15.50 


-15.62 


-15.72 


0.1 


2 


-12.02 


-13.16 


-17.93 


-18.08 


-19.36 


-20.02 


0.1 


3 


-14.65 


-15.66 


-17.95 


-18.01 


-18.88 


-19.50 


0.1 


4 


-17.61 


-18.04 


-19.10 


-19.10 


-19.66 


-19.91 


0.2 


1 


-16.00 


-15.88 


-15.81 


-15.67 


-15.69 


-15.85 



-15.73 -15.32 -16.65 -17.86 -19.01 -19.10 -20.83 -21.62 

-20.78 -20.52 -21.75 -22.50 -22.67 -22.65 -22.77 -22.64 

-20.20 -19.94 -21.22 -22.02 -22.23 -22.22 -22.80 -22.93 

-20.22 -19.82 -20.65 -21.12 -21.28 -21.27 -21.67 -21.68 

-15.91 -15.38 -16.59 -17.70 -18.89 -18.98 -20.80 -21.59 

-20.70 -20.44 -21.80 -22.49 -22.86 -22.89 -23.01 -23.02 

-20.15 -19.87 -21.23 -22.01 -22.45 -22.46 -22.98 -23.11 

-20.28 -19.87 -20.80 -21.24 -21.50 -21.52 -21.91 -21.91 

-16.18 -15.49 -16.55 -17.58 -18.76 -18.86 -20.75 -21.58 



22.49 


-23.53 


-27.67 


-17.96 


-18.40 


-19.70 


-20.60 




-22.57 


-22.59 


-22.50 


-20.06 


-19.32 


-17.42 


-15.87 




-24.57 


-26.06 


-29.56 


-20.04 


-19.69 


-22.04 


-22.60 




23.13 


-24.57 


-28.31 


-18.95 


-18.42 


-20.43 


-21.47 




-22.44 


-23.52 


-27.39 


-17.92 


-18.37 


-19.64 


-20.42 


38.23 


22.86 


-22.89 


-22.82 


-20.28 


-19.69 


-17.74 


-16.19 


38.23 


23.72 


-26.16 


-29.41 


-20.20 


-19.82 


-22.09 


-22.48 


38.23 


-22.37 


-24.69 


-28.29 


-19.15 


-18.61 


-20.47 


-21.38 


38.23 


22.39 


-23.52 


-27.30 


-17.87 


-18.35 


-19.61 


-20.49 


39.87 



Note. — The electronic table supplies the absolute magnitude of each SED template (AGN: 1, E: 2, Sbc: 3 and Im: 4) for a broad range of photometric bands as a function 
of redshift, along with the distance modulus DM. A complete version of this table can be found in the electronic edition of the journal. Templates have been normalize d to a 
"bolometric" luminosity of 10^^ Lq. The absolute magnitude we present here corresponds to the canonical definition of the absolute magnitude (as in, for example, e qn. 2 6|Of lHogel 
Il999l ) plus the K correction term. This allows the calculation of photometric redshifts and K corrections from the table. This is equivalent to Tables 3 and 4 of I AOS! (j^e their 
caption for details). All magnitudes are in the VEGA system except for NUV, FUV, z and the WISE bands, which are in the AB system. Note that no template inc(9^orates 
IGM absorption and the AGN template is unreddened. For calculating photometric redshifts and /C-corrections allowing for variable IGM absorption and obscuration of tile AGN, 
public numerical routines are provided on-line at http://www.astronomy.ohio-state.edu/~rjassef/lrt. 
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Table 3. Photometric Redshift Estimation Summary. 



Sample 


(yzl{\^z) 


Az 


68.3% 


95.5% 


99.7% 


Median 


All 


0.155 


0.057 


0.049 


0.240 


0.721 


-0.020 


Point Source AGNs 


0.261 


0.184 


0.174 


0.544 


0.777 


-0.121 


Extended AGNs 


0.149 


0.050 


0.049 


0.175 


0.812 


-0.006 


Galaxies 


0.128 


0.041 


0.042 


0.127 


0.429 


-0.017 



Note. — Summary of the photometric redshift calculations for our main 
catalog and each sub-sample discussed in §4.21 The table shows for each case 
the value of cj^/ (l+z) (as defined by eqn. [1]), Az (the 95% clipped distribution 
C2/(l+-2)), the ranges of — 2;s|/(l-|-z) encompassing 68.3%, 95.5% and 99.7% 
of the distribution and the median value of Zp — Zg. 
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0.1 1 10 

A(/2m) 

Fig. 1. — {Upper panel) The wavelength coverage of our data set [black), as described in 
§21 The sohd hne shows the coverage when photometry upper hmits are included while the 
dashed line only considers detections. We also show both histograms for the sample used to 
construct the galaxy SED templates alone {gray) which excludes AGNs (see §3.ip . {Bottom 
panel) Filter sensitivity curves for all bands included in our data set are shown by the solid 
black lines. Filter sensitivity curves for the upcoming WISE mission are shown by the dark 
gray solid lines at half-height to distinguish them from the other bands. 
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Fig. 2. — Edwards representation of a four set Venn diagram showing the number of AGNs 
selected by each photometric criteria and the overlap between them. The top panel shows 
this for our full sample while in the bottom panel we have restricted it to objects with z > 1 
in order to show how the photometric selection criteria overlap in a contamination free, but 
limited, sample. Note that in each panel there are four geometrical figures, each representing 
a photometric selection criteria: MIPS (magenta rectangle). X-ray (blue rectangle), radio 
(green circle) and IRAC (red "peanut" shaped region). The intersection between the regions 
show the number of objects targeted simultaneously by each combination of criteria. 
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Fig. 3. — Resulting best fit templates from applying the algorithm of §3.11 to the data 
set described in ^ (b lack solid) compared to their initial guesses {black dotted) and to the 
templates derived in lAOSl {gray solid). Note that the bottom axis of each panel shows 
wavelength while the top shows frequency. The bars at the bottom of each panel show the 



rest-frame wavelength coverage of our data set {black, see Fig. [T]) and of that used in |A08 
{gray). 
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Fig. 4. — Converged AGN SED template using the iRichards et al.l (l2006[ l template as a 
starting point {solid line) and our standard mod el {dotted line). The Figure also shows the 
original SED template of IRichards et al.l ( 120061 ) {dashed line). The templates are poorly 
constrained near > 20/im due to the lack of powerful low redshift quasars. Note that the 
bottom axis shows wavelength while the top shows frequency. 



- 34 - 




Fig. 5. — Average LAGN/-Z^Host obtained using the converged iRichards et al.l ( 120061 ) AGN 
SED template as a function of the ratio found using our standard AGN template. On the 
solid line the ratios are equal, while on the dashed lines they differ by a factor of 10. Error 
bars show the dispersion in each bin of LAGN/-^Host determined with the standard AGN 
template. 
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Fig. 6. — Best fit reddening to the AGNs in our sample as a function of the contribution of the 
AGN component to the total bolometric luminosity. The AGN luminosity has been corrected 
for the estimated extinction. The top and middle panels show the objects targeted as AGN 
candidates by AGES according to their photometric properties (see §2]). The bottom left 
panel shows the reddening distribution for objects classified as AGNs by the spectroscopic 
pipeline while th e botto m right panel shows it for the objects classified as Type 2 AGNs by 
Moustakas et al.l (120091 ) . In each panel, the solid black line shows the median E(B-V) as a 
function of LAGN/-^^Host while the dotted and dashed lines show the 85% and 95% contours 
respectively. 
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Fig. 7. — Best fit SEDs for four typical AGNs {top four panels) and for the AGN with the 
90**^ percentile worst in the sample {bottom panel). The solid line shows the best fit SED, 
which is composed of a non-negative combination of our optimized AGN {dotted)., E {dashed), 
Sbc {long-dashed) and Im {dot-dashed) templates. In each panel we indicate the reddening 
apphed to the AGN template. The flux of each filter, when detected, is shown by the sohd 
triangles, while the open squares show the expectation from the best-fit. Arrows mark upper 
limits when they are available. The AGNs shown in the top four panels have SEDs typical 
of, from top to bottom, a Type 1, a Type 2, a very low luminosity AGN and a Type 1 with a 
very strong host component. For each object, the of the fit and the spectroscopic redshift 
measured by AGES is shown on the upper right corner of the corresponding panel. 
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Fig. 8. — Best fit SEDs for four typical galaxies (to'p four panels) and for the galaxy with the 
90*^^ worse xl the sample {bottom panel). The different line styles and point types follow 
the same convention as in Figure [71 In the top four panels we show galaxies dominated by, 
from top to bottom, the elliptical, Sbc and Im templates, and a combination of all three. 
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iX 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

Spectroscopic Redshift 

Fig. 9. — Comparison of photometric and spectroscopic redshifts for all the galaxies in our 
sample (right) and for the 95% with the most accurate redshift estimates (left). For a given 
value of the photometric redshift, we show the median spectroscopic redshift of a galaxy with 
that Zp estimate (dotted line), and the 68.3% (solid line) and 95.4% (dashed line) range of 
the spectroscopic redshift distribution. We only show Zp values for which there were enough 
objects to determine the 68.3% and 95.4% hmits on both sides of the median. 
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Fig. 10. — Comparison of photometric and spectroscopic redshifts for the extended ito'p) 
and point source AGNs {bottom) in our sample considering all (right) and the 95% with the 
most accurate redshift estimates (left). Different line styles have the same definition as in 
Figure O The open circles show individual objects with either Zg > 4 ot Zp > 4. 
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Fig. 11. — All best fit SEDs (gray) and their average (black) for point source AGNs with good 
photometric redshifts {\zp — Zs\ < 0.2) in four ranges of the AGN reddening. The number 
of objects used in each case is shown in the top-right corner of each panel. All SEDs are 
normalized to match at 3.6;Lim. The vertical lines show the wavelength of Lyman-a (dotted), 
4000A (dot-dashed) and 3.6/xm (dashed). The extreme scatter at wavelengths shorter than 
Lya is due to large variations in IGM absorption with redshift. 
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Fig. 12. — Same as Figure [TTl but for point source AGNs with bad photometric redshifts 
(|zp — Zs| > 0.5). Note how these sources have very flat, featureless SEDs compared to those 
in Fig. [m For most of these AGN, the Lyman break is not in the range of the data. 
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Fig. 13. — Difference between the photometric and spectroscopic redshift of all point source 
AGNs in our sample as a function of the bolometric luminosity ratio between the best fit 
AGN and host components. Arrows show upper or lower limits when either the AGN or the 
host components were not used in the best fit. When the AGN is more luminous than its 
host {rightward of vertical gray dashed line), the photometric redshift accuracy is poor. 
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Fig. 14. — IRAC color-color diagram of all point source AGNs in our sample. Black points to 



the objects with bad photometric redshifts {\zp 



> 0.5) while the rest are shown in gray. 



Large points mark objects classified by the modified SDSS pipeline to have active nuclei, 
while the rest correspond to the small symbols. Notice that the majority of solid circles lie 



at [3.6] - [4.5 



of Stern et al. 



colors redder than 0.6. The solid lines mark the AGN color selection region 
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Fig. 15. — Comparison of the bolometric luminosity ratios between the AGN and host com- 
ponents for all point source AGNs {solid gray points). Points with black borders correspond 
to objects with bad photometric redshift determinations {\zp — Zg] > 0.5). When either the 
AGN or the host component is not detected, an arrow shows the upper/lower limit of the 
respective ratio. Note that the ratios are not systematically biased when the photometric 
redshifts are used instead of the spectroscopic ones. 
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Fig. 16. — IRAC band color-color diagram of / < 21 .5 SDWFS sources [gray dots). Overlaid 
are the AGN selection region of IStern et al.l (120051 ) {solid black boundaries) and the color 
tracks (as a function of redshift) of our E {dot-dashed black line), Sbc {dashed black line) 
and Im {solid black line) galaxy SED templates, and of our AGN template for no reddening 
{solid gray line) and for a reddening of E{B — V) = OA {dashed gray line). For the galaxy 
templates, color tracks are shown for redshifts between 0-3, while for the AGN template 
they are shown between 0-10. Each heavy dot marks an increase of unity in redshift for the 
galaxies and an increase of 2 units of redshift for the AGN. For each template, the heavy dot 
surrounded by a circle marks z = while those surrounded by a square mark the terminal 
redshift {z = 3 for galaxies and z = 10 for the AGN). 
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Fig. 17. — Same as Figure ITBlbut in the colors used by lLacy et al.l ( 12004J ) to define their AGN 
sel ection crite r ion (s olid black boundaries). We also shown the updated selection criteria 
of iLacy et al.l ( 120071 ). modified to limit the problems with contamination by low redshift 
galaxies {dotted black selection boundaries). Line and point styles have the same definition 
as in Figure [161 
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Fig. 18. — WISE band color-color diagram of all SDWFS sources that would be detected in 
the [3.3] and [4.6] bands {small dots), or in all four channels simultaneously {large dots). Each 
color represents objects that were targeted for spectroscopy by AGES (but not necessarily 
observed) as IRAC {red), MIPS {magenta). X-ray {blue) and/or radio {green) AGNs, or as 
non-active galaxies {gray). These are applied sequentially, so the MIPS points are only those 
that were not IRAC selected and so forth. The selection boundaries using all four channels 
{solid gray line) and just the [3.3] — [4.6] color {dot-dashed gray line) are also shown. For 
reference, we show the color tracks of the SED of Arp 220 {dotted black line) and of our E 
template {dot- dashed black line) as a function of redshift. For each of the SEDs color tracks, 
heavy dots mark an increase in redshift of unity, while the large open circle marks z = and 
the large open square marks z = 3 for Arp 220 and z = 2 for the E template. 
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Fig. 19. — Color tracks of our SED templates and of the SED of Arp 220, overlaid on top 
of all SDWFS sources that would be detected at [3.3] and [4.6] by WISE {gray dots). The 
different line-styles have the same definition as in Figures fT6] and ITSl Note that, because of 
the shallow depths of the WISE mission, we only show the color tracks up to 5; = 2, 6 and 
3 for galaxies AGN and Arp 220 respectively. Most of the Arp 220 track extends off the red 
edge of the Figure (see Fig. [18] for the complete color track). 
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Fig. 20. — This Figure shows the effective i-band magnitude Umits of the WISE survey 
for unreddened AGNs (top) and moderately reddened AGNs {E{B — V) = .1 , bottom) , for 
pure AGN SED on the left and the average AGN + Host SEDs of IRAC-selected AGNs 
on the right. Each black line shows the i-band magnitude corresponding to the depth of 
each WISE band as a function of redshift. The bold sohd (bold-dashed) black line shows 
the i-band magnitude limit when detection is required in all four (two shortest wavelength) 
WISE bands. The dotted gray lines show the SDSS hmits for z < 2.3 {i < 19.1) and z < 2.3 
{i < 20.2). 



